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ABSTRACT 

We continue our earlier studies of quasi-periodic oscillations (QPOs) in the 
power spectra of accreting, rapidly-rotating black holes that originate from the 
geometric "light echoes" of X-ray flares occurring within the black hole ergo- 
sphere. Our present work extends our previous treatment to three-dimensional 
photon emission and orbits to allow for arbitrary latitudes in the positions of the 
distant observers and the X-ray sources in place of the mainly equatorial posi- 
tions and photon orbits of the earlier consideration. Following the trajectories of 
a large number of photons we calculate the response functions of a given geome- 
try and use them to produce model light curves which we subsequently analyze 
to compute their power spectra and autocorrelation functions. In the case of an 
optically-thin environment, relevant to advection-dominated accretion flows, we 
consistently find QPOs at frequencies of order of ~ kHz for stellar-mass black 
hole candidates while order of ~ mHz for typical active galactic nuclei (~ 1O 7 M ) 
for a wide range of viewing angles (30° to 80°) from X-ray sources predominantly 
concentrated toward the equator within the ergosphere. As in our previous treat- 
ment, here too, the QPO signal is produced by the frame- dragging of the photons 
by the rapidly-rotating black hole, which results in photon "bunches" separated 
by constant time-lags, the result of multiple photon orbits around the hole. Our 
model predicts for various source/observer configurations the robust presence of 
a new class of QPOs, which is inevitably generic to curved spacetime structure 
in rotating black hole systems. 

Subject headings: accretion, accretion disks - - black hole physics - - X-rays: 
galaxies — stars: oscillations 
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1. Introduction 



One of the more interesting findings of the X-ray timing analyses of the light curves of 
accreting compact objects in the past two decades has been the discovery of quasi-periodic 
oscillations (QPOs). These are broad features in the power spectra of these sources that 



range up to kHz frequencies in neutron star low-mass X-ray binaries (see, e.g., Ivan der Klis 



2000) and up to ten s a nd even hundreds o f Hz in accreting black hole binary systems (see, 



e.g.. lCui et all 1 19981 and Istrohmaveill2001 al ibi), including QPO pairs with 2:3 frequency com- 
mensurability in the power spectra of certain black hole candidates (e.g. XTE J1550-564 
and GRO J1655-40). Besides galactic binary sources, the presence of QPOs has also been 
reported in at least on e ultra- luminous X-ray sources (ULX), namely NGC 5408 X-l (see 
Strohmaver et al. 2007f) at v ~ 20 mHz; in the context of active galactic nuclei (AGNs) 
Gierlihski et al.l (120081 ) recently discovered a ~ 1 hour X-ray periodicity in the narrow-line 
Seyfert 1 (NLS1) RE J1034+396, in which the periodicity is apparent in its light curve in 
distinction with those of the galactic binary black hole systems mentioned above that do not 
exhibit any such behavior. 

The nature of QPO phenomenon in accreting black hole systems is at present largely 
unknown; it is also possible that this phenomenon is not the manifestation of a unique pro- 
cess, but that the observed diversity of QPO frequencies and sites may be due to a multitude 
of processes specific to the particular frequency and site. For the case of QPOs associated 
with black hole systems, the absence of an underlying, rotating solid surface object, implies 
that plausible explanations involve by necessity processes associated with the surrounding 
accretion disks. These include, among others, accreti on disk precession due to the frame^ 
dragging effects of a ra pidly-rotating black hole (e.g 



Schnittman et al. 20061 ) (also see, e.g., Stella k, Vietri 



Qui et al. 1998; Merloni et al. 1999 



19981 for the same effect invoked to 



acco unt for the kHz QPOs of Low Mass X-ray Binaries), accretion disk oscillat ory modes 
(e.g. iNowak et al.lll997l : lAbramowicz fc Kluzniakll200ll : lKatoll200ll : iDonmea 120071 for disko- 
seismology), or the distribution of X-ray emitting "blobs" spanning a limited range aroun d 
some specific radius of the Keplerian disk surrounding the black hole (e.g. iKarad 119991 ). 
There have also been studies of the correlations of the QPO fre quencies with the properties 
of the associated X-ray spectra (e.g. iTitarchuk fc Fioritol l2004j ) . 



An altogether different notion to QPO origin has been put forward recently by lFukumura fc Kazanas 
These authors proposed that high frequency QPO (HFQPO) could 



( 120081 . hereafter FK08 

be produced as a consequence of light echoes of X-ray flares occurring within the ergosphere of 
a Kerr black hole; they showed that due to frame- dragging of individual flare photons, a finite 
number of them reach far away observers after an integer number of additional orbits around 
the black hole to produce a geometr y-induced "light echo" of the original flare (see, also, 



Meyer et al.l 120061 ; iBursa et al.l 120071 . for a similar discussion) . Studying two-dimensional 



(2D) near-equatorial photon propagation in a fully general relativistic calculation, FK08 
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showed that dragging of inertial frames leads approximately 15% of the photons to a specific 
direction at infinity (an observer) after a time lag of At ~ 14 M (where M is black hole 
mass), independent of the relative position between the observer and the photon source. 
The independence of the lag on the relative phase between the observer and the source 
then guarantees a second peak in the Autocorrelation Function (ACF) at lag r ~ 14M, 
which according to Fourier analysis, results in a prominent HFQPO signal at a frequency 
z/qpo — 1/At ~ 1.4 (IOMq/M) kHz (where M is the solar mass), even for a light curve that 
consists of flares randomly distributed within the ergosphere. Since the closest an accretion 
disk can approach the black hole is the radius of Innermost Stable Circular Orbit (ISCO), 
this effect is possible only when the ISCO lies within the black hole ergosphere, a condition 
that limits this effect to black holes with dimensionless spin parameter a/M > 0.94. 

The study of FK08 showed the prominent role of frame- dragging in producing lags 
independent of the relative source - observer position; as shown there, the absence of frame 
dragging for flares in a Schwarzschild geometry leads to lags which depend on the relative 
observer - source position and, for random flare positions, to absence of QPO features in the 
corresponding power spectra. 

The original study of FK08 outlined the fundamental aspects behind this work, namely 
the importance of frame-dragging and the persistence of QPO presence, even for random 
positions of the X-ray flares within the ergosphere. However, it was restricted in its scope 
in that it was constrained to observers and sources near or on the black hole equatorial 
plane, i.e. 6" s ~ ^ ~ 7r/2. Realistic constraints demand an enlarged study to include, at 
a minimum, observers at larger latitude positions, not least because it is expected that the 
column density of sources for lines of sight near the equator are sufficiently high to preclude 
propagation of photons without interaction with the surrounding matter; this could then 
lead either to their absorption or to their scattering and introduction of additional, random 
lags to their trajectories erasing the otherwise potential QPO signal. The goal of the present 
work is to remedy this deficiency by expanding the study of the corresponding lags to three- 
dimensional (3D) geometries allowing thus for arbitrary latitudes of both the observers as 
well as of the source of photons. 

We would like to stress that the present models of "geometry-induced" QPOs (and also 
those of FK08) are no t meant as an account of the so-c alled HFQPOs presented in the 
literature to date (e.g. Cui et al. 1998 ; Strohmayer 2001al Jbl) . For one thing, the predicted 



QPO frequencies are 5 — 10 times higher than the observed frequencies mentioned above 
(for the same black hole mass), and, for another, they do not produce their 2:3 frequency 
commensurability. Instead, the predicted QPOs should be viewed as an altogether new 
feature indicating the intrusion of an accreting gas within a black hole ergosphere, thereby 
providing an unequivocal evidence for the presence of rapidly rotating black holes. 

Our paper is organized as follows: In §2 we provide a general description of our 3D 
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QPO model, the details of the photon kinematics and response function of the system, 
and a prescription for constructing stochastic model light curves. In §3 we compute the 
autocorrelation functions (ACFs) and power spectral densities (PSDs) of the model light 
curves to show that they generally exhibit the QPO features as anticipated. Finally, in §4 
we review our results, make contact with observations along with the limitations of our model, 
and conclude by discussing prospects of future work. Finally, in the Appendices we discuss 
the details the notion of locally isotropic emission in Kerr geometry and provide explicit 
functional forms relating the photon impact parameters to their local emission angles. 



2. Model Details 

The present treatment follows closely that of FK08, but generalizing the photon tra- 
jectories to 3D, thus allowing for positions of both the source and the observer at latitudes 
other than near equatorial. So, we consider photon emission in Kerr geometry of dimension- 
less angular momentum of a black hole a, described in Boyer-Lindquist coordinates (r, 9, (ft). 
As in FK08, here too we employ the usual geometrized units (G = c = 1 where G is the 
gravitational constant and c is the speed of light) and normalize distance and time with the 
black hole mass M. For example, 1M corresponds to 1.5 x 1O 5 (M/M ) cm in distance and 
5 x 1O _6 (M/M ) sec in time. We consider the observer position at an arbitrary location 
(r ,9 ,4> ), while correspondingly denoting the source coordinates as (r s ,9 s ,(p s ). This ar- 
rangement, therefore, allows for source locations not only at random azimuths but also at 
random polar angles. In thin accretion disks, the flares are likely confined to small latitudes, 
i.e. 9 S ~ 7r/2, but f or quasi-spherical ac creting flows such as advection-dominated accretion 
flows (ADAFs; e.g. iNaravan fc Yilll994h . the overall source latitude can be much higher. As 



discussed in following subsections, in this last case, one can, in addition, weigh the emission 
at each latitude to fit one's preferred flow model. 



2.1. X-ray Flares 



The X-ray emission of the sources we consider is regarded as the incoherent sum of 
a large number of instantaneous (i.e. short compared to the dynamical time) flares. We 
assume the typical duration of such a flare to be ~ 1M, presumably originating from local 
magnetic reconnections (e .g. ISvensson fc Zdziarskil Il994f) or standing s hocks in relativistic 



magnetized accretion (e.g. iTakahashi et al.ll2002l; iFukumura etal.l 120071) . Concerning gener- 



ation of X-ray flares within the ergosphere, iKaras fc Kopaceki (120081) have at tributed it to 
magnetic reconnection induced by frame-dragging, while iKoide fc Arail (120081 ) have consid- 
ered magnetic reconnection within the ergosphere as a means of extracting energy from the 
black hole rotation. In general, these X-ray flares do not necessarily have to be confined 
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Fig. 1. — (a) Schematic view of the system considered in this paper, (b) Poloidal projection 
of the geometrical setup in this problem. 



exactly to the equatorial plane, and for this reason, in this work, we relax this restriction 
on the source vertical location. Therefore, in the specific cases that the emitting plasma is 
optically thin, the observed signal will contain contributions from sources in both upper and 
lower hemispheres. 

It is important to mention that in this model both polar and azimuthal positions (9 S , <p s ) 
of each X-ray flash are randomly assumed (as long as they lie within or close to the ergo- 
sphere) while their cylindrical radius r s sin8 s is fixed for simplicity to the radius of the ISCO 
of a Kerr hole with a/M ~ 0.99, i.e. to r s sin# s ~ 1.455M. One of our primary interests 
is thus to probe the dependence on the source and the observer latitudes as illustrated in 
Figure [Tb. To avoid the introduction of artificial coherent signals associated with orbital 
motion of these emitting sources, we also consider the flares to be randomly distributed in 
time in the following fashion: each X-ray burst occurs with a Poisson distribution of mean 
value equal to the Keplerian orbital timescale; AT = T orb | ln{rnd(0, 1)}| where T orb is the 
orbital time of the source at (r s , 6 S ) and rnd(0,l) is a random number between and 1. We 
define that T orb = 2ii{{r s sin# s ) 3 / 2 + aM 1 / 2 }/!/ 1 / 2 . 

The orbiting X-ray emitting plasma is assumed to be in Keplerian motion, either below 
or above an equatorial plane at 9 S ^ where % denotes the i-th source, i.e. rotating with the 
local Keplerian angular velocity (modified to take into account the black hole spin a) 



s [ (r sin 0) 3 /2 + a MV2 
with all the quantities being evaluated at the source position. In the lab frame or the locally 



(1) 
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non-rotating reference frame (LNRF) each source has a total three- velocity of 

~Asin9 



SA 1 / 2 



(2) 



where we define A = r 2 — 2Mr + a 2 , £ = r 2 + a 2 cos 2 0, A = (r 2 + a 2 ) 2 — a 2 A sin 2 9, and 
frame-dragging of the local inertial frame is denoted by uj = IMraj A. In the LNRF, as one 
can see, rotation of a local inertial frame u has been subtracted by definition. 



2.2. Null Geodesies 



Neglecting external interactions between photons and particles, each null geodesic (i.e. 
photon orbit) is uniquely characterized by two constants of motion, £ and rj, where £ is 
the axial component of angular momentum a nd rj closely related to its polar component 
(e.g. iBardeen et al.lll972t IChandrasekharl Il983l ) . For given rj and £, a photo n orbit in Kerr 



geometry is generally gove rned by the following equations of motion (e.g. 
19721 ; IChandrasekhaj[l983h : 



Bardeen et al. 



r 2 £ 2 
9T? 



(r 2 + a 2 )E + aA£ - a 2 A sin 2 9 
AS 

H 2 - A{ v + (£ - a) 2 } , 
V + - £ 2 esc 2 6) cos 2 9 , 
a(2Mr - ag) + £A esc 2 6 



(3) 

(4) 
(5) 

(6) 



where S = r 2 + a 2 — a£ and an overdot denotes derivative with respect to the affine parameter. 
Note that equations (j3J) and (jSJ) provide only the squares of r and 9. In order to avoid the 
issue of determination of appropriate signs for r and 9, for the actual numerical determination 
of the orbits we use instead the second derivatives of equations (TjJ and (jSJ) 



9 



2rfZ - r (r - M) {q + (£ - a) 2 } - E |2rr - a 2 sin(2#)^ | 

rS 2 

£ 2 cot 9 esc 2 ^ - a 2 sin(29)9/2 - S |2rr - a 2 sin(2fl)^} ^ 2 

^E 2 



(7) 



The initial conditions necessary for the integration of the orbit equations involve, besides 
the source position, (r s , 9 S ), also the two angles of photon directions and the values of r and 
9; the last four quantities are not independent but they can be expressed in terms of (r s , 9 S ) 
and the values of the impact parameters 77 and £ as discussed in §2.3. 
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Fig. 2. — Geometry of local emission in the LNRF described in the text. 



2.3. Locally Isotropic Emission 



In the absence of any information concerning the angular distribution of the photons 
emitted at a given flare, the most conservative assumption is that the emission is isotropic 
in the rotating fluid rest frame. However, the proximity of the X-ray sources to the black 
hole and the associated source motion (v s ~ 0.6 c) and strong field geometry along with the 
fact that the coordinate system for the integration of the photon orbit equations is that of 
the LNRF, require that the notion of photon isotropy be transformed to this frame and that 
the impact of the Kerr metric on the notion of isotropy be considered with some care. Below 
we discuss our prescription of local isotropic emission, while more details can be found in 
Appendix A. 

The geometry of photon emission is depicted in Figure [2j We assume the black hole 
spin axis to be in the z-direction while the local fluid velocity v s in the y-direction; ip is the 
angle between the emitted photon and the instantaneous direction of source motion (y-axis) 
in the fluid frame, while ip' is the corresponding angle in the LNRF. For a given ip the precise 
direction of emission is determined by the angle \ confined in the plane orthogonal to the 
y-axis. The relation between ip and ip' is then (the angle \ remains invariant) 



where (3 = v s /c, while the differential solid angle and opening angle are transformed respec- 
tively as 




(9) 





(10) 



Isotropy in the fluid frame implies equal number of photons per unit solid angle N 0l i.e. 
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Fig. 3. — Comparison between the fluid rest-frame and the LNRF with [3 = 0.5 for (a) 
Differential photon number distribution per emitting angle as a function of respective emit- 
ting angle in the fluid frame (dark) and the LNRF (gray), (b) Corresponding bin-corrected 
photon distribution (i.e. weighting factor) in the fluid frame (dark) and the LNRF (gray). 
We take JV = 1000 and Aip = tt/50. 

dN/dVt = No over the entire sky in the fluid rest frame. However, when computing photon 
geodesies we choose incrementally not the solid angle but the angles ip and \- Then, the 
number of photons emitted per increment of the polar angle ip is dN/ dip = 2nNo sin ip and 
therefore the increment of the angle \ should be equal to l/(27riVnsin'0) so that each solid 
angle element at a given ip, x contains one photon. 

The Lorentz boost to the LNRF skews the photon distribution with angle resulting 
in larger number of photons being emitted along the direction of source motion with the 
distribution as a function of the polar angle ip ( or given by the function 

tfo /a _ (! -/5 2 )sin^' (l + /3cos^)sin^ _ 

W>*)= (l-^eos^) 2 = (1-^ = fiM) ■ (U) 

In Figure [3^l the photon distribution as a function of the polar angle ip f° r P — (dark dots 
for the fluid frame) and also as a function of ip' for [3 = 0.5 (gray dots for the LNRF) is 
shown. We take iVn = 1000 and Aip = 7r/50. However, as indicated above in equation ffTUj) . 
the Lorentz boost transforms also the width of the the angular bin A^; because of that, 
as shown in Appendix A, the number of photons per angular bin is invariant (as it should 
for photon number conservation), proportional to sin-^A^. So each angular bin Aifj' in 
the LNRF has the same number of photons as the corresponding bin A^ in the fluid rest 
frame, but these photons are now received at a smaller angle ip', given in terms of ip by 
equation (J9j). This distribution is shown in Figure [3b, where the number of photons per 
angular bin f(f3, ip')Aip' = (dN/dip')Aip' = (dN/dip)Aip - for the same values of N Q and Aip 
as in Figure [3^, - is shown as a function of the angle ip (dark dots for the fluid frame) and 
ip' (gray dots for the LNRF). It is apparent that a given number of photons per bin is found 



- 9- 



at smaller angles in the LNRF rather than in the fluid rest frame. 

For our photon orbit calculations we transform an equally spaced (discrete) set of angles 
ipi in the source rest frame to the corresponding angles %p[ in the LNRF and release at each 
such bin a number of photons equal to 2ir No sin ipi. These photons are then distributed each 
in a single of as many equal intervals in the x angle. The next step involves expressing local 
emitting angles 5' and 7' (both measured in the LNRF) used in the integration of the photon 
trajectories in terms of the angles xp' and \. The relations between these sets of angles can 
be obtained from trigonometric relations of the geometry shown in Figure [2] as 

cosip' = sin 5' sin 7' , cos 5' = simp' cos x , (12) 

where the first of these equations was obtained from a spherical triangle for the cosine of 
the angle ip' while the second one from relations of the projection of the geometry of the 
same figure onto the plane perpendicular to the y-axis (i.e. in the shaded ellipse of the same 
figure). 

The next step in the calculation takes int o account the effects of the local geometry 



curvature in the definition of local isotropy (see iFukumura &: Kazanas 1 120071 . as an example 
of this issue in the simpler case of the Schwarzschild geometry) by relating the angles 5' 
and 7' to the photon four-velocity components. In the LNRF the photon's local emitting 
polar and azimuthal angles (5', 7') are related to the four-velocity components through the 
relations 

A 1 / 2 r 

v r = — -. — r = sin 8' cos 7' , (13) 
At 

A\ 1/2 9 

-) ,=00.*, (14) 
1 Asm 9 , 

v* = ir - u\ = sin 5' sin V , (15) 



EAV2 I i 



' A \ 1/2 f) 

cos8> = [-) -j , (16) 



which lead to 



a; t 

C ° W = AVWHmBi-ui' (17) 

evaluated at the emitting location (r s ,9 s ,0). Using equations (jSJ) through (jfJJ) we can elim- 
inate t, f, 9 and <p from the above equations ( |T6l) and ( jTTl) in favor of the two impact 
parameters 7] and £. We are thus led to two relations that connect the photon emission 
angles (5', 7') to the photon impact parameters (£,??). In fact, one can invert these relations 
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Fig. 4. — (a) Locally isotropic photon trajectories from a source at (sin# s , <p s ) = (89°, 0°) for 
an opening angle of ip = 80° (tp' = 46°). We set a/M = 0.99 and N = 100. (b) Observed 
(cumulative) photon distribution N tot (9) / N tot (90°) as a function of covering angle measured 
from the pole for various source geometries. It is normalized by the total photon number 
over tc/2. We assume an equatorial source of 9 S = 89° (dark curve) and mid-latitude source 
of 53° (gray curve), in comparison with stationary flat-space case (dotted curve). 

to obtain (£, rj) as a function of these angles, i.e. £ = £(S',j') and rj = tj(5',j') as shown in 
Appendix B. This now completes the problem of the choice of the orbit parameters consis- 
tent with isotropic emission in the fluid frame: For a given set of equally spaced values of ipi 
(locally isotropic source in the fluid rest frame) one can compute and also determine the 
number of photons to be emitted and values of the azimuthal angle bins \i\ these are then 
used to compute the angles 5' and 7', and then the corresponding impact parameters rj and 
£ along with the initial values of f(5', 7') and 9(5', 7') that are used in the integration of the 
photon orbit by equations ([7j) and (jSJ). 

An example of 3D photon trajectories from a locally isotropic emission is shown in 
Figure Hk. where we consider emission with a half-opening angle of if) — 80° (ip' = 46°) 
from a source at (9 s ,(p s ) = (89°, 0°). Throughout this paper we assume a/M = 0.99 as in 
Figure UK, unless otherwise stated. Recall that our model considers now X-ray sources from 
both the upper and lower hemispheres with respect to the equator because of the optically- 
thin assumption of surrounding matter (see Fig. [1]). Hence, the net observed emission will 
consist of a sum of these two components. We have confirmed that these photon trajectories 
are indeed almost symmetric with respect to the equator as expected. It should be noted, 
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however, that for some particular emission angles this symmetry in trajectory with respect to 
the equator can be largely broken even for a small deviation of the source position from the 
equator. Due to the high velocity of the source (/3 ~ 0.6) for the values of r s and a chosen, 
most of the emission is beamed forward in the direction of instantaneous source velocity. In 
fact we find that many backward-emitted photons (i.e. of ip > 90°) subject to the Doppler 
and frame-dragging effects of the fast black hole rotation are forced to corotate in the hole 
direction (see FK08 for details in 2D cases). We also find that backward-emitted photons 
with t/> > 140° end up crossing the the horizon, making no contribution to the observed 
signal. We should mention that the photon trajectories considered here also include many 
multiple orbits and higher-order photons; i.e. we also collect those photons that orbit around 
a black hole multiple times with some time-lags. The physical significance of these time-lags 
will be further discussed in §2.4. 

To illustrate the effects of Doppler beaming and the black hole's strong gravity for 
the considered source location (r^sin^ ~ 1.45 M), we present in Figure Hb the cumulative, 
normalized distribution of photons iV tot (^)/iV tot (90 o ) received by a distant observer at r/M = 
100, as a function of the observer's inclination angle o . In the figure the solid dark curve 
corresponds to a nearly equatorial source {6 S = 89°), the gray curve to a source at the edge 
of the ergosphere (8 S ~ 53°) for the fixed cylindrical radius of the source, while the dotted 
curve to the flat-space distribution. It is apparent in this figure that photons emitted from 
equatorial position (in particular) are directed primarily at low latitude directions. While in 
flat-space 50% of the photons (horizontal dotted line) are emitted between angles 0° — 60° 
as expected, in the equatorial source case of 6 S = 89° the half-way point has been extended 
to 9 ~ 82°. Similarly, while in flat-space polar angles 0° — 45° include ~ 30% of the entire 
photons emitted, in the present Kerr geometry they include roughly 10 times less. This is 
an important issue related to the equivalent width of the (broad) fluorescence Fe lines whose 
red wing is presumably emitted by plasma closest to the ISCO of fast-rotating black holes 
with spin parameter similar to that used in producing our Figure Hb. 



2.4. Response Functions and Light Curves 

Following a methodology similar to that discussed in FK08, we first compute a source's 
response function seen by a distant observer whose radial and azimuthal positions are set 
in this work to (r Q /M, <p ) = (100,0). This is a reasonable choice because at this distance 
the photon orbits are already nearly straight and thus their final (polar and azimuthal) 
directions are almost the same as those in flat spacetime. For a given source position at 
(@s,i, 4>s,i) around a rapidly-rotating black hole we collect isotropically emitted photons (in 
the local rest-frame of an orbiting source) keeping track of its arrival time t and final position 
(8, <p) for each photon. As a result we can construct a response function seen by the observer 
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Fig. 5. — Synthetic response functions from a single X-ray source at 9 S = 89° with (a) 
4> s ~ 4° and (b) ~ 180° seen by an observer of 9 = 60°. Arrows indicate a constant lag of 
~ 1AM. Values of the other parameters are the same as in Figure HI 



Ii{t, 9 Q ; 9 St i, 4> s ,i), which tells us photon counts as a function of photon's arrival time from the 
i-th X-ray source at (9 S ^, <fi s ,i)- For an optically-thin environment we have contributions from 
sources at different source latitudes for the same source azimuth, however, for simplicity we 
assume that the sources closer to the equatorial plane are dominant, although it is found 
that this assumption does not make much difference to our final results shown §3 below. 

The temporal resolution in the response function is taken to be 1M (corresponding to 
5 x 10~ 5 sec for a 1OM black hole), therefore the signal received by the observer is uniformly 
binned with this numerical accuracy. For computational purposes we approximate this signal, 
as precisely as possible, by a set of narrow Gaussians, also of width 1M. Note the response 
function from each single source contains a characteristic temporal profile due to distinct 
photon trajectories between the source and the observer, and the profile also strongly depends 
on the observer's polar position 9 Q as well as the source position (9 S ^, 4> s ,i)- The observer's 
collecting angle is chosen to be A9 D = ±1°. To produce a statistically significant outcome 
even for observer's polar angles 9 Q <, 30° we consider a large number of sampling (a few 
10 6 ) per source. However, the common feature they share is found to be a constant time 
separation of ~ 14 M between the major peaks of individual response functions, as found in 
FK08. To clearly illustrate these features of the response functions, we exhibit in Figure [5] 
the response functions from flares at (a) <ft s ~ 4° and (b) ~ 180° for 9 S = 89° measured by an 
observer of 9 = 60°. A constant time lag of ~ 14M (equivalent to ~ 0.7 msec for a 10M Q ) 
is clearly present. Due to multiple orbits around a fast-rotating black hole some photons 
make extra (integer-multiples) full orbits before reaching the observer. 

The response functions are different in the precise positions of the major peaks, as 
one would expect, but, as argued above, the time-lag between the major peaks is roughly 
constant, At ~ 14 M, regardless of the source location (9 S , <p s ). The reason for this behavior 
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was discussed in detail in FK08 (it remains the same in the present case too): it is due to 
the frame-dragging effects which force all photons emitted at any relative position between 
source and the observer, to reach the latter only by "swinging" around the black hole; in 
following such trajectories, a non-negligible fraction (~ 15 — 20%) of them reach the observer 
after multiple orbits around the black hole, inducing "light echoes" of nearly constant lag as 
shown in Figur e El This is in agreement with similar relativistic calculations of broad iron 



line profiles by IBeckwith fa Pond (120051 ) where they find that these higher-order photons 
contribute maximally ~ 20 — 60% of the total luminosity of the system. The persistent 
presence of such a constant time-lag in the response functions is the key factor leading to 
QPO features in the power spectra shown in §3. It should be noted that in the absence of 
strong frame- dragging (i.e. sources outside the ergosphere or Schwarzschild instead of Kerr 
black hole), the time interval among the peaks in the response functions is phase- dependent 
and precludes the presence of QPO in the power spectra for random source positions (see 
FK08). 

As discussed in §2.1 the X-ray interval between flares AT is Poisson-distributed with 
mean value comparable to the Keplerian period at the ISCO. Therefore, the collective re- 
sponse function from the ensemble of X-ray flares (i.e. the bolometric light curve) is obtained 
by summing over every contribution from all the sources at {0 s ^,<j) St i). That is, 

n 

i(t,e o ) = J2 I i(-t,8o;0s, i ,<f>s > i) , (is) 
i=i 

where n is a total number of random X-ray sources in our simulations. As an example, a 
synthetic light curve is shown in Figure [6^ (i.e. superposition of all the response functions 
similar to Fig. [5]) and its first 20Mx msec segment in Figure [6b where Mi = M/(10M & ). 
In this specific case, the number of flares considered is n = 2000 with the source positions 
confined to near the equatorial plane (8 S = 89°) while the observer is set at o = 60°. Because 
of the random positions (4> s ^) and bursting times of individual flares, the model light curve 
appears to be (in fact it is constructed to be) totally random with no apparent coherence 
(even in the absence of background noise). Nonetheless, the characteristic time-lag among 
the peaks in flux implicit in the form of the response function in Figure [5] is encoded in the 
light curve and it can be extracted (if present) by applying the standard time-series analysis 
techniques discussed in the next section §3. 



3. Power Spectra and Autocorrelation Functions 

In this section we discuss in detail the timing analysis metrics, ACFs and PSDs, of the 
model light curves produced using the prescriptions discussed above, pertaining to rapidly- 
rotating black holes; more specifically we study the dependence of the ACFs and PSDs on 
the sources' position and geometric structure as well as the position of the observer. 
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Fig. 6. — (a) Model light curves constructed from the synthetic response functions similar 
to those in Figure We assume n = 2000 random X-ray sources from 9 S = 89° for 9 Q = 60°. 
(b) The first 20Mi msec segment of (a). Values of the other parameters are the same as in 
Figure 0J 

As shown in §2 above, the form of the response functions and hence of the model light 
curves we produce depend mainly on the sources' latitudes (given that we average over 
their azimuths <p s ,i)- Because the effects of frame-dragging are most prominent for sources 
near the equator, we first examine cases where most of the X-ray emission is concentrated 
preferentially near the equatorial plane, i.e. of 9 Sj i = 89° (and corresponding positions in 
the lower hemisphere). Regarding the observer's inclination angle o we consider a wide 
range from nearly face-on (9 D ~ 30°) positions to nearly edge-on (9 ~ 85°) ones. Since 
our primary goal in this investigation is to generalize the previous 2D results of FK08, we 
restrict ourselves to a fast-rotating black hole case with a/M = 0.99 for which the source 
X-ray photons are subject to strong dragging of inertial frame (i.e. inside the ergosphere). 
Note that for a/M = 0.99 the ISCO radius is naturally well inside the ergosphere near the 
equatorial plane. 

We also consider, for a given observer's angle 9 Q , the dependence of the ACF and PSD 
on the latitudinal distribution of the sources within the ergosphere. To this end we compute 
also the response function for source positions at different vertical heights with A9 S = 9° by 
dividing the polar angle into five representative zones; i.e. 9 S = 89°, 80°, 71°, 62°, and 53° (see 
Fig. [lb). We then examine the effects of the sources' vertical extent by computing light curves 
from sources with incrementally larger vertical extent, i.e. from sources whose emission 
consists of the sum of point sources at a number of heights corresponding to the following set 
of inclination angles: 9 S = 9$9° (i.e. equatorial sources) , 9s9<=+so°, #89 o +80 o +7i°, #89 o +80°+7i°+62°, 
and 6 , 89°+8o°+7i°+62°+53°- We do not consider sources at heights that would correspond to 
angles smaller than 53°, because for the given value of black hole spin (a/M = 0.99) and 
a vertical source whose equatorial foot point is at the ISCO, a larger height would place 
that section of the source outside the ergosphere, thereby producing no contribution to the 
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QPO features via frame-dragging effect. Given that the effects of frame dragging decrease 
with increasing source height we expect the QPO effects to decrease with increasing vertical 
source extent. 

The timing properties of the model light curves are explored for different relative posi- 
tions between the source's and observer's polar angles by computing the ACF, R(t;9 ,9 s ), 
and the corresponding PSD, P{y\ 9 Q , 9 S ). For a discrete light curve I(t , 9 D \ 9 S ) the normalized 
ACF is given by 

p/v.fl a\ - I(t k ,9 ;9 s )I(t k+j ,9 ;9 s ) 

where L is the number of time bins in the light curve of I. The corresponding PSD, in units 
of (rms/mean) 2 /Hz, is given by 

2T 

P{yf % e m e.) = w \a \\ (20) 

ph 

where T is the time duration of the light curve I and aj is the Fourier amplitude defined by 

1 L 

I(t k ,9 ;9 s ) = -J2^ m/L > (21) 

i=i 

and iVph is the total photon counts in the entire light curve of /. Also, Tj and Vj are timescale 
and frequency of a characteristic QPO signal, respectively. 



3.1. Dependence on Viewing Angle 

In this subsection we first present in Figure [7] a series of ACFs (left panels) and PSDs 
(right panels) as a function of the observer's inclination (or viewing) angle 9 Q for model light 
curves from both near-equatorial sources, i.e. of 9 S = 89°, and a source that is vertically 
extended in height, modeled as a sum of individual sources at heights that correspond to 
polar angles 9 S = 89°, 80°, and 71°. These are shown in Figures [Th^ and b (equatorial source) 
and c and d (vertically extended source). 

As noted in FK08 the QPO signal in the light curves is imprinted by the frame dragging 
effects induced by the black hole spin. It is seen that a characteristic timescale of r ~ 
14M = 0.73Mi msec is persistently present for a wide range of inclination 9 , which is a 
manifestation of the constant time-lag in the response function in Figure [5j Since these 
effects are more prominent for photon orbits that lie near the equator, this signature is more 
pronounced for arrangements that maximize the propagation of the received photons in this 
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Fig. 7. — Autocorrelation functions (left panels) and power spectral densities (right panels) 
for various observer's inclination angles where 9 S = 9 8 q° (upper panels) and 6 S = 89 o +8QO+71 o 
(lower panels). We show 6 = 30°, 40°, 50°, 60°, 70°, 80°, 85° from bottom to top curves with 
Av = 0.6 Hz bin size. Vertical positions of the curves are arbitrarily shifted for presentation 
purpose. Values of the other parameters are the same as in Figure HI 



region. As shown in Figures [7^ and b this effect is the most optimized for o ~ 50° — 60°; 
for observers at smaller inclination (polar) angles the effect is weaker because the influence 
of frame-dragging is much reduced for photons that propagate close to the black hole spin 
axis. This becomes apparent by the decrease of the ACF peak at r ~ 0.73M! msec, or 
equivalently, the corresponding QPO amplitude with the decrease in O (Fig. [7H and b). The 
QPO amplitude decreases also for 8 Q > 60°; the reason for that is the increase of the fraction 
of the poloidal plane orbits that connect the source and the observer; these photons do not 
orbit around the black hole with coherence and as such they do not contribute to the QPO 
amplitude but in fact they dilute it. That is, the total signal is the result of a superposition 
of a mixture of both coherent and incoherent signals. 

Similarly, as shown in Figures and d, a vertically extended source leads to reduction in 
the corresponding amplitudes of the ACF and QPO since the frame-dragging effects decrease 
as the source approaches the edge of the ergosphere. Comparison of Figures [7h and b with 
Figures [Tfc and d, which exhibit respectively the ACFs and PSDs for a set of observers at 
inclinations 9 Q = 30°, 40°, 50°, 60°, 70°, 80°, 85° for an equatorial (a,b) and a vertically 
extended source (c,d), confirms these notions; the QPOs are more prominent for sources 
confined near the equator than for those which are vertically extended. The former also 
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persists over a wider range of values of o , being strongest for O ~ 50° — 60°, while the 
latter, being intrinsically weaker, they disappear for o > 60°. 

Finally, while the QPO and ACF secondary peak amplitudes do depend on the observer 
inclination, the corresponding QPO frequency and ACF peaks remain constant respectively 
equal to vqpo ~ 1.4M-f kHz and r ~ 0.73M! msec, a value representing the length of the 
photon paths around the ISCO, becoming considerably wider only for 6 Q < 30°. 



3.2. Dependence on Source Concentration 

If the X-ray sources originate from processes other than those associated with the surface 
of a thin disk (e.g. flaring events due to magnetic reconnection, polar shocks or wind/jets), 
then the source may not be exclusively distributed in an equatorial fashion. In this subsection 
we examine in greater depth the effects of extending the source dimension vertically on the 
ACFs and QPOs respectively. To this end we have produced the light curves starting with 
an equatorially concentrated source (9 S = 89°) and consider the light curves from more 
extended sources by incorporating into the light curve also photons from sources at larger 
heights, or smaller values of the angle 6 S ; this is done incrementally by considering additional 
sources of the same intensity at the following set of source polar angles: 8 9°+so°, #89 o +80 o +7i°> 
6*89°+80°+7i°+62°, and 6 l 89° + 8o o +7i o +62°+53° • We do not consider sources at any larger heights 
because for the conditions used here they would be outside the ergosphere and they would 
not contribute to the QPO signal by frame-dragging. 

In Figures and b we present respectively the ACFs and PSDs for an observer at 
9 = 30° as the source height progressively increases to 6 S = 53° (i.e. sources extending all 
the way up to the static limit). Overall, the QPO signal does not appear to be very strong 
but it is present. We see that the additional (incoherent) photons dilute those with constant 
lag to cause the QPO to disappear with the source height. Similarly, in Figures Eb and d 
we present the ACFs and PSDs for an observer with high inclination of 6 = 60°. We see 
that in this case the QPO becomes again prominent when the X-ray sources are distributed 
near the equator as the photons that reach these directions are under a strong influence of 
frame-dragging. It appears that the latitudinal extension of the sources up to 6 S < 71° can 
yield a discernible QPO signature, although it is relatively weak. 



4. Summary and Discussion 

In the present work we have extended the earlier 2D treatment of FK08 to consider the 
3D geometry of photon emission on the timing properties of accretion flows in the vicinity 
of rapidly spinning black holes. This has allowed us to explore the effects of placing both 
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Fig. 8. — Autocorrelation functions (left panels) and power spectral densities (right panels) 
for various source distributions where o = 9 30 o (upper panels) and o = 6o ° (lower panels). 
We show 6 S = # 89 o, 6 S go +80Oj 6 S9 o +80 o +71 o, 6 , 8 9o +80 o +71 o +62 o, and #89°+80°+7i°+62°+53° from bot- 
tom to top curves with Av = 0.6 Hz bin size. Vertical positions of the curves are arbitrarily 
shifted for presentation purpose. Values of the other parameters are the same as in Figure HI 



observers and sources at latitudes much different from the equatorial ones considered in FK08 
and, using these results, to consider also photon sources which are extended, i.e. sources 
whose vertical dimension is comparable or larger than their radial one. Our results are 
fundamentally consistent with those of FK08: we find that for the same reasons discussed in 
FK08, here too, frame-dragging causes the photons to follow trajectories that by and large 
lead to a significant fraction of them reaching the observer after an additional orbit around 
the black hole; as such, for a single X-ray flare, an observer sees multiple bunches of photons 
that arrive with a constant time-lag of order of ~ 14M (i.e. light echoes), independent 
of the source position relative to the observer. It is the independence of the lag on the 
source position that leads to the QPO in the signal of even a random ensemble of X-ray 
bursts, provided that they take place within the ergosphere. As elaborated in FK08, the 
resulting QPOs are different in character to those examined so far in the literature in that 
they do not require any modulation/oscillation as such in the signal. The difference in their 
character can be easily assessed in the data by calculating, in addition to the PSDs, also the 
ACFs, which in the present case has the form shown in Figures [7] and [HI i.e. a peak near 
zero (self-coherence) with a second one at r ~ 1AM = 0.73Mi msec (QPO), very different 
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from that of an oscillatory form, corresponding to QPOs due to some oscillation in the X- 
ray flux (Fukumura, Dong, Kazanas, & Shrader, in preparation). In summary, for sources 
concentrated near the equator (e.g. 6 S <, 85°), we find the presence of persistent QPOs (and 
a series of harmonics) of vqpo ~ 1-4 kHz for a 1OM black hole (equivalently ~ 1.4 mHz for 
a 1O 7 M AGN) for a wide range of line-of-sight viewing angles (from ~ 30° to ~ 80°). This 
QPO signal is the direct consequence of photon trajectories undergoing multiple rotations 
around the black hole due to its strong frame-dragging effects near the equator. As the 
source extends toward the mid-latitudinal region (~ 53°) the QPO signal is weakened by 
more dominant incoherent signals. 

The QPOs proposed in this work should be viewed as a new class of QPOs in view of their 
(i) expected frequency bandpass and (ii) their underlying mechanism, generic to the dragging 
of inertial frames on individual photons. As such, we do not believe these are associated 
with the observed HFQPOs often detected in accreting black hole systems (whose frequency 
is lower by a factor of ~ 5 and some of which exhibit a 2:3 frequency commensurability 
that the present model does not provide for). The fundamental requirement for the presence 
of QPOs of the type discussed herein is that the X-ray flares would take place within the 
ergosphere of a rapidly-rotating black hole. The requirement that the ISCO lies within the 
ergosphere, then, implies that the black hole spin be a/M > 0.94. While this appears to be 
a rather strong constraint, one should point out that broad Fe line fits by emission from thin 
accretion disks around Kerr black holes, provide con sistently estimates of the spin parameter 
near a/M ~ 0.99 (e.g. iBrenneman fe Reynoldsll2006l ). This value is consistent with that used 
in the present work. 

The above results are relevant if the accreting gas surrounding the central engine is 
relatively tenuous and optically transparent to radiation from these X-ray sources. For 
example, this should be the case for optically-thin ADAFs where the accreting plasma dis- 
tribution is quasi-spherical (though most of emission is still originating from the midplane 
of the gas), in contrast to the equatorial structur e of a standard, thin accretion disk (e.g. 
Shakura fe Sunyaevlll973l : iNovikov fe Thorndl 19731 ). Such geometrically thick, optically thin 
configur ations appear to be associated with objects such as our Galactic Center source Sgr 



A* 



(e.g. 



Yuan et al.ll2003l ; iMeyer et al.ll2006l ). and Low Ionization Nuclear Emission-line Re- 



gion sources or LINERs (e.g. NGC 4258; see Quataert 1999) and nearby giant ellipticals 
in low-luminosity AGNs or LLAGNs (e.g. M87; see Di Matteo et al. 2003). If the central 
emission region in these sources is indeed optically-thin characterized by ADAFs or its more 
generalized class of radiatively- inefficient accretion flows (RIAFs; see, e.g. Yuan et al. 2003), 
then one may expect a transparent environment in which photons can escape the production 
region without much additional scattering. 

QPO features have recently been associated with the PSD of AGN: XMM-Newton ob- 
servations of the the bright nearby Seyfert 1 galaxy Mrk 76 6, known to be a NLS1 , exhibited 
statistically significant peak in the PSD of its light curve (IMarkowitz et all 120071 ). modeled 
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as Lorentzian of vqpo ~ 0.46 mHz, similar to QPOs detected from other black hole binary 
systems. Assuming an estimated mass of the Mrk 76 6 nucleus (~ 5 x 10 6 M Q ), suggested 
from optical H/3 emission line velocity dispersion (e.g. I Wandell 120021 ) . we can estimate the 
frame-dragging QPO frequency to be vqpo ~ 2.7 mHz, i.e. a factor of ~ 5 higher than the 
observed frequency. It is possible that other NLSls with slightly more massive nuclei (order 
of 10 7 M Q ) may exhibit a QPO frequency that can be explained in the context of our model. 
For example, a radio-loud NLS1, PKS 0558 -504, shows a rap id X-ray flare presumably with 



a heavier nucleus of M ~ 4.5 x 1O 7 M (e.g. IWang et al.ll200lh wh i ch ma y make this source a 



good candidate for testing our prediction. Recently, iMeyer et al.l (120061 ) discussed a charac- 
teristic time-lag due to higher-order emission from our Galaxy Center (Sgr A*), which may 
also be a promising target for our investigation (because of larger mass and longer timescales 
compared to stellar-mass black holes). In addition to Seyfert Is above, perhaps intermediate 
mass black holes possibly a ssociated with ULXs ma y also be promising sites to search for the 
QPOs proposed herein (e.g. lStrohmayer et al.ll2007l ). Although it may not be observationally 
easy to disentangle a potential QPO signal from high luminosity continuum (presumably at 
near-Eddington rate) from these objects, it remains to be studied. 

Our discussion so far has been confined to X-ray flare sources that orbit the black hole 
in the azimuthal direction with no significant radial motion. One could argue, however, that 
any source in the vicinity of an accreting black hole is likely to be plunging-in radially at a 
good fraction of c. This is certainly the case for sources interior to the ISCO. As pointed out 
above, in order that a source be both at the ISCO and within the black hole ergosphere, the 
spin parameter a should be a/M > 0.94. For more moderate values of the spin parameter, 
say, a/M ~ 0.9, the ISCO is outside the ergosphere and plasma traversing the latter should 
have considerable radial velocity in addition to its azimuthal motion. Since there is no 
reason for which this plasma could not produce X-ray flares (for example, shocked plasma in 
plunging regions can in principle serve as a good candidate for such a thermal/nonthermal 
X-ray sources, see Fukumura et al. 2007), the entire analysis described in this work can be 
extended to this circumstance too. It is expected, however, that because of its (substantial) 
radial velocity, the emitted photons would be Doppler beamed in the radial direction too and 
a large number of them would be lost through the black hole horizon; as a result, the QPO 
production efficiency may decrease to an unobservable level for sufficiently small values of 
the spin parameter a. A detailed study of this arrangement and the limiting value for which 
this approach can lead to QPO signals is deferred to future work. 

A related issue is that of the detectability of QPOs produced in the way described 
above. In this work we neglected the background distribution of signal/noise associated 
with other physical processes (e.g. emission from accreting flows and/or corona), photon 
counting statistics (Poisson noise) and also instrumental responses to the signal. Depending 
on the QPO signal strength relative to the externally contaminating (e.g. accreting gas 
and coronal) emission intensity, the induced light echo can be sufficiently weakened to a 
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statistically indiscernible level. To assess a potential observability of the predicted QPO 
signal in this scenario we need to combine the synthetic signal with the stochastic noise. If 
fraction of the observed X-rays does not participate in the production of the light echoes 
discussed herein, one should add to our synthetic signal an underlying constant (stable) 
component with some signal-to-noise ratio. Finally, in order to gauge the observability of 
the signal we propose, one should also add the noise associated with photon statistics and 
detector background. Such a detail treatment of the observability of the QPOs proposed in 
this work is deferred to another future work. 

We would like to thank the anonymous referee for a number of useful and insightful 
suggestions. K.F. is grateful to Fotis Gavriil for helpful discussions about timing analysis. 
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5. Appendix A 

In this Appendix we describe in details the notion of local isotropy and derive the 
equations necessary to calculate the subsequent photon trajectories in Kerr geometry. 

In order to correctly describe local isotropy in emission from a source orbiting at rela- 
tivistic speed, we first consider the differential photon distribution as a function of the angle 
ip between the source direction and the photon emission in the fluid frame 

dN dN d{cosip) 
dip d(cosip) dip 

The notion of local isotropic emission is the requirement that the number of photons per 
local solid angle should be the same for all the direction of emission, i.e. that 

where N denotes the number of photons within the solid angle dfl which must be conserved 
from one frame to another and dQ = 2k sinipdip. Therefore, we obtain 

dN 

— = 2nN smiP, (24) 
dip 

which provides the weighting factor as a function of ip. 

In the LNRF (or lab frame), where the photon trajectories are computed, we similarly 



find 

dN dN d(cosip) d(cosip') 
dip' d(cosip) d(cosip') dip' 

where again 

dN 



(25) 



= 2ttA . (26) 
d(cosip) 



This preprint was prepared with the AAS IATgX macros v5.2. 
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Using equation (J9J), we can express equation (j25|) as 

dN 



27rJV o /0M): (27) 



where 



_ (1+^cosj0) sim/> 
W = (i_ £2)i/2 ' 

Note that dN/dip' -> dN/dip as /3 -> and we can see J* {dN/ dip') dip' = /* (dN/dip)dip = 
47iNo regardless of the value of /? as expected. Hence, f(j3,ip) describes the differential 
photon distribution per unit opening angle (i.e. the weighting factor) in general inertial 
frames. 

One can easily convert the above distribution into the differential photon distribution 
over a finite bin Aip' (in the LNRF) given by equation fllOp . Note that Aip' is not constant 
and one finds 

M ' = {%) ^ = 2nN ° ^ ' (29) 

which states that the number of photons emitted within the corresponding angular bin in 
the two frames is indeed the same (i.e. independent of inertial frame), as one expects. Then, 
it is guaranteed that for the entire sky 

» = t O ^; = t (i) a * - ^ [ *"»* = ^« . (») 

where the integer i represents a discretised grid point in the polar opening angle and n 
denotes the total number of bins of the same angle from to n. This prescription of photon 
distribution allows one to compute the number of photons in each local opening angle ip', 
consistent with isotropic emission at the fluid frame. For each such polar angle, we then 
launch photons in the azimuthal (about the source velocity) direction x (see Fig. [2]) between 
and 27r in equal intervals of width 

AX = 2n (jhp' A ^ ) ~ iV sin ipAip ' ^ 

with this prescription, then, guaranteeing the proper number and photon directions in the 
LNRF consistent with isotropic emission in the fluid rest frame. 

Besides the Lorentz boost of the angles to the LNRF one must in addition relate these 
angles to the photon angular momenta (£,77), as these are the quantities needed in the 
integration of the photon geodesies. A simple geometry of this problem relates the opening 



-25 - 



angle ip' and the azimuthal angle x t° t ne corresponding polar and azimuthal angles, 5' and 
7', as described in §2.3 

cos ip' = sin 5' sin 7' and cos 5' = sin ■?// cos x > (32) 

which is schematically illustrated in Figure[2l 



6. Appendix B 



In Appendix A we discuss how to prescribe local isotropy from an arbitrary X-ray source 
position in Kerr geometry, which enables us to calculate the emission direction (5', 7') of an 
individual photon. In order to subsequently integrate the geodesic equations (photon trajec- 
tories) per photon given by equations ([3]) through (J5J) one needs to translate such a directional 
information into the photon's angular momentum (equivalent to impact parameters) in the 
LNRF. Since we are dealing with 3D null geodesies each photon is characterized by two 
constants £ and rj where £ is the axial component of angular momentum and rj is closely 
relat ed to the polar component of angular mo mentum pg by 77 = p 2 e + £ 2 cot 2 9 — a 2 cos 2 9 



(e.g. iBardeen et al.lll972l ; IChandrasekharl Il983l ) . Here we derive analytic expressions in the 
most general case for (£,77) in terms of the local emission angles (5', 7'). 

Considering equations (1IB]) and ( 1TT1) where the left-hand-side is related to the local 
emission angles and the right-hand-side to geometry and velocity, one can eliminate the 
derivatives (t, r, 9, <fi) using equations (j5J) through to express the right-hand-side in terms 
of (£, 77). By solving these equations one can obtain analytic forms of £ = £(5', 7') and 
rj = r](5', 7') explicitly as 

1 - /j 2 VAC sin 2 5'± {(1 -/i 2 )AS 2 C 2 sin 2 6' esc 2 7 '} 1/2 " 



^CA{£ 2 



AS 2 esc 2 i - (1 - ,u 2 V 2 A 2 sin 2 5' 
- a 2 (l - /1 2 )} + (1 - /i 2 ) (C - loAO 2 cos 2 5' 



(33) 



(34) 



;i-^ 2 )ca 

where ( = r 4 + a 2 {/i 2 A + r(r + 2M)} and /1 = cos 6 1 with all the quantities to be evaluated at 
the source position (r s , 9 S ,4> S ). The plus/minus sign for the second-term in the numerator of 
equation ( |33l) depends on the initial direction of photons emitted in such a way that £ > 
for all emission; i.e., plus for photons initially emitted (locally) in the forward direction 
(0 < 7' < 7r) and minus for backward direction (vr < 7' < 2tt) with respect to the LNRF 
(which can be in rotation relative to the observer). Note that there is an "offset" factor in 
directional information due to this rotation of local inertial frame [due to the quantity u in 
eqns. (13"3"|) and (1541) ]. In the absence of frame-dragging (a/M — > 0) we see that 

( 1 — /Lt 2 \ 1/2 

£ — > ±r ( - - - t ) sin 5' sin 7' , (35) 



1 - 2M/r 
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which is positive for both forward and backward emission with sin 5' > 0, as expected (see 
Fig. [2]). This also reproduces correct axial component £ in Newtonian geometry (r — ► oo); 
when emission is exactly radial, either away (7' = 0) or toward (7' = it) the rotation 
(symmetry) axis, £ will vanish as required. When emission is completely parallel to the 
rotational axis (either 5' = or it), one also recovers £ — > 0. 

As a special case we can see that for emission from an on-axis (z-axis) source, i.e. as 



indicating the absence of axial angular momentum, as anticipated. For another special case 
where emission is confined in the equatorial plane or (x,y)-plane as 9 — > tt/2 (/i —>■ 0) and 
5 f — > 7r/2 we see 



e -»• (/w 1) 



£ — > and r] — > finite , 



(36) 



£ — > finite and 77 — > , 



(37) 



where 77, which in this case is exactly p% 



<g, vanishes, also as expected. 



